Local-energy density functionals for an N-dimensional 
electronic system in a magnetic field 
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We present a general approach for the construction of 
the exact local-energy density functionals for a uniform N- 
dimensional electronic system in a magnetic field. For arbi- 
trary dimension, we obtain explicit expressions for the matter, 
kinetic, and exchange density functionals. In the zero-field 
limit, we recover the usual N-dimensional Thomas- Fermi the- 
ory. As an application of our results, we develop a current- 
density-functional theory, in the spirit of the Thomas-Fermi- 
Dirac approximation, for an inhomogeneous many-electron 
system in a magnetic field. 



I. INTRODUCTION 

It is now well established that a density-functional ap- 
proach to the ground state and dynamical properties 
of many-electron systems is a viable alternative to the 
usual wavefunction approach involving the solution of the 
single-particle Schrodinger equation. The proof by Vig- 
nale and RasoltE (VR) that the scalar and vector poten- 
tials characterizing a many-electron system are uniquely 
determined by the single-particle electron and current 
densities (the so-called current-density-functional theory 
CDFT) has placed density-functional theory on a rigor- 
ous footing for the description of many-electron systems 
involving electric and magnetic fields. 

In this work, we present a general approach for the 
construction of the local-energy density functionals for an 
N-dimensional electronic system in a homogeneous mag- 
netic field. Unlike the Kohn-Sham single-particle scheme, 
in which only the exchange-correlation functional is un- 
known, this statistical formalism requires the knowledge 
of both the kinetic and exchange energy functionals in 
the presence of electric and magnetic fields. Our ap- 
proach is to consider the exact canonical density ma- 
trix for a noninteracting N-dimensional electron gas in 
a homogeneous magnetic field. The single-particle den- 
sity matrix is obtained via the inverse Laplace transform 
of the canonical density matrix, as originally formulated 
by March and Murrayo. The diagonal component of the 
single-particle density matrix yields the exact density dis- 
tribution of the system. The exact kinetic and exchange 
energy-density functionals may also be explicitly eval- 
uated from the single-particle density matrix. Although 
these functionals are obtained for a homogeneous system, 
an application of the local-density approximation (LDA) , 
in the spirit of Thomas- Fermi-Dirac theoryB (TFD), al- 
lows these functionals to be also used in inhomogeneous 



systems. 

Our paper is organized as follows. In Sec. [n|, we 
define the N-dimensional canonical density matrix for a 
noninteracting electron gas in a homogeneous magnetic 
field, and discuss some of the mathematical properties 
of the associated N-dimensional single-particle density 
matrix. Then, in Sec. [II, we consider the evaluation 



of the kinetic and exchange energy-density functionals 
for arbitrary dimension and magnetic field strength. Ex- 
plicit forms for the kinetic and exchange energy-densities 
are obtained for three-and-two-dimensional systems. In 
Sec. [V, we present an outline for the construction of a 



current-density-functional theory based on the TFD ap- 
proximation in two-and-three-dimcnsions ; making use of 
the density functionals obtained in Sec. III. Finally, in 
Sec. g we offer our concluding remarks. 



II. DENSITY MATRIX OF A N-DIMENSIONAL 
HOMOGENEOUS ELECTRON GAS IN A 
UNIFORM MAGNETIC FIELD 

Long ago, Sondheimer and WilsonB considered the 
evaluation of the diamagnetism of free electrons. The 
central tool in their calculation was the determination of 
the canonical density matrix, which in 3-dimensions (3D) 
is given by (including spin degeneracy) 



C(r,r';/3) = 
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where the symmetric gaugeB has been used and uj — 
eB/mc. For purely formal reasons, we extend the defini- 
tion of C(r, r'; (3) to N-dimensions in the following way 
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with (N = 2,3,...). The single-particle density matrix 
in N-dimensions is then obtained via the usual inverse 
Laplace transform relational 



P(r, r') 



e-f-ioc 



df3 e 



0s f 



2ni 13 



-C(r,r';/3) , 



(3) 



with Ef being the Fermi energy, and the contour chosen 
to be closed in the left-half of the /3-plane (see Fig. [l]). 
We defer an exact treatment of Eq. (^|) until the next 
section. 





FIG. 1. (a) The contour to be used when d is integral, and 
(b) when d is half-integral. 

The electron number density p (spin-averaged, zero 
temperature) is immediately given by taking the trace 
of the density matrix in Eq. (|3j), viz., 



, M(7) . 2(s) v« (7 >, 



with 



F [d \l) 



= 2 7 s 



2wi s sinh(s) 



(4) 



(5) 



and for notational convenience, we have defined s = 
?iuj{3/2, 7 = £f/(huj), and d = N/2. Depending on the 
dimension of the space, d can be either integral or half- 
integral. 

In the case that d is half- integral (i.e., odd-N), the in- 
verse Laplace transform defined by J 7 ^' always contains 
a branch point in the complex s-plane at s = 0. There- 
fore, care must be taken when evaluating J 7 ^ (7) at the 
branch point. The contour chosen for the evaluation of 
Eq. (||) is shown in Fig. [l]b, where a cut has been in- 
troduced along the negative real axis, and the large arc 
does not pass through any of the poles of the integrand. 
Following the method of Sondheimer and Wilson p we 
find that the integral in Eq. (|J) can be written as a Rie- 
mann zeta function plus the sum of the residues at the 
poles Sfc = ±kni 



ds 



=2 7 s 



2-ni s d sinh(s) 
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s^sinh(s/c), 



, (6) 



where a denotes the contour FEDC in Fig. |l]b. The first 
integral in Eq. (|^) can be separated into two parts by 
using an asymptotic expansion of the integrand at the 
branch point 
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where 
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k\ ds k \sinh(s) / s=0 
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The first integral is evaluated by means of Hankel's for- 
mula for the gamma function, viz., til 



ds 2ls _ z (27)*- 1 

- — e ' s — 

2m 



T(z) 



(9) 



and the second integral can be cast as real since the part 
contributed by the circle around the origin (see Fig. |l]b) 
is zero. The sum of the residues is readily computed, and 
we finally have 



^ (d) (7) = £ 
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where 
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The first term in Eq. ( [L0[ ) gives the TF expression for the 
number density in a homogeneous magnetic field. The 
second and third terms are purely quantum-mechanicaL 
and are the origin of the de Haas- van Alphen effect. 
These oscillatory terms are clearly a result of the exact 
treatment of the canonical density matrix, aaid are not 
obtained within the usual TF approximation!] 

Integral values of d allow for an exact treatment of Eq. 
(||). Using the contour shown in Fig. [j]a, T^ d \^f) can be 
expressed in the form 
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= 2 7 s fc 



s^sinh(sfe) 
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where s k — kiri (k — 0,±1,±2, ...). The k = contri- 
bution represents the usual TF approximation, whereas 



2 



the terms associated with the poles at ±kwi are respon- 
sible for the quantum-mechanical de Haas-van Alphen 
effect. As in the half- integral case, 

jr(d)( 7 ) for even-N 
will be the sum of steady terms (s& = 0) plus oscillatory 
terms (sk — ikni): 



2(-l) fe 
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where the terms are given by 

k=0 d\ ds d Vsinh(s) ) 



(13) 
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In the limit to — + 0, the asymptotic behaviour of 
^f( d ) (Vy) f or arbitrary d is given bytJ 



(27)° 



(15) 



and r 
viz. J? 



r(d + i)' 

reduces to the N-dimcnsional TF theory. 



(d) _ 
Ptf — 
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For large magnetic fields, the general behaviour of 
j:( d ) f or half- integral d is difficult to obtain using the 
present treatment. Once again, the cause of this diffi- 
culty lies in the branch point at s — of the complex 
s-plane. We can circumvent this problem by writing Eq. 
(0) in the following formO 



ds e (27s-2fc-l)s 

2m s d 



where we have used the series expansion 
1 



sinh(s) 



k=0 



-(2fc+l)s 



(17) 



(18) 



for the sinh(s) function. Only those terms satisfying 2k + 
1 < 2 j survive the integration in Eq. (|l7j), and since 
max [27] = 2, we have 



^)( 7 ) = T ^(2 7 -1) 



d-l 



(19) 



Equations (Q) and (|T^) provide the high-magnetic field 
electronic number density for arbitrary values of d. 



given by Eq. (J3J) . However, in the present context, a sim- 
plification for the evaluation of the kinetic energy density 
is possible. Our first step is to notice that the Heavyside 
function may be represented as the-inverse Laplace trans- 
form of the exponential function£3 



e+ioo 



dj3 e^ £ f-^ 
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The kinetic energy density £ki n may then be expressed as 

e kin = (r I H0{e f -H)\t) 

- lco d(3 (r I He^ £ f-^ | r) 



2m 
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For arbitrary dimension, e^in takes the form 
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where 



27 



(d+l)/d 



(22) 
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Equation ( |22j ) is valid for arbitrary magnetic field 
strength. 

In the limit of zero magnetic field, the asymptotic be- 
haviour of H^("f) for arbitrary d can be deduced from 
Eq. @, viz., 

hWm r(d + 2)[r(d + i)]^-[r(d + i)]^ 

H (7) = WT2) ■ (24) 



Thus, using the above value of H ^ (7) in Eq. , we re- 
cover theJcinetic energy density of the N-dimensional TF 
theory.B£3 As the dimension of the space goes to infinity, 
we obtain 

ff M (7) = l, (25) 
whereby the kinetic energy density becomes 



III. ENERGY-DENSITY FUNCTIONALS 

The kinetic and exchange energy functionals can both 
be calculated from the single-particle density matrix, 



Substituting Eq. ( |l~9| ) into Eq. 
limit of # (d) (7)> viz., 



(26) 

gives the high-field 
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# (d) (7) 



4 7 r(d)^r(rf + i)(2 7 -i) c 
r(d + 1)2^(27- 1)^ 
2(27- i) rf r(d)^ 



r(d + i)2^(2 7 -i)- 

In the limit d — > 00, we obtain 

#W( 7 ) = 2 7 , 
and the kinetic energy behaves as 



(27) 



(28) 



A. N=3 



The electron number density in three-dimensions can 
be obtained from Eqs. (Q) and ([To|): 



0(3/2; 
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(27) 



3/2 



1 



^ (3/2) (7)+f:^^cos(2/,,- 



7T 



6v27T7 

T 



. (34) 
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The other density functional of interest is the N- 
dimensional exchange energy e cx . According to the usual 
prescription, e cx is given byu 



' ' / - / I P(r,r') 



"4 ' * 



(30) 



where the integration is understood to be over the Tri- 
dimensional volume. The exchange energy density has 
a much richer s-plane topology than the kinetic energy 
density because of the essential singularities arising from 
the coth(s) term in Eq. (|2|) . These essential singularities 
make it difficult (if not impossible) to obtain a useful, 
closed analytical form for the N-dimensional exchange 
energy for arbitrary magnetic fields. Nevertheless, some 
analytical progress can be made if we make use of the 
following expansions 



exp(— acoth(s)) 
sinh(s) 



2^L„(2a)e 



-a -(2n+l)s 



(31) 



where L n (a) = Ln\&) is the Laguerre polynomial. It is 
then easy to show that the N-dimensional single-particle 
density matrix has the form 



P(r,r') 



'-(x'y-y 



' x >]TL„(2a)e 
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2(d-l) 

muj ^ 

5 = 7F 2^^ 



z[f . (33) 



Equation (|32| ) in Eq. (|30|) yields, in principle, the N- 
dimensional exchange energy for arbitrary magnetic field. 
Of course, the comments made in Sec. |l| with regard to 
the electron number density also apply to Eq. (|32]). 

Below, we illustrate the utility of our formalism 
through specific examples in the physically relevant 3D 
and 2D systems. 



where Q^^^i^j) is given by Eq. (|Tl]). In the zero field 
limit, we can make use of Eq. dl^) to get 



,(3/2) 
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3tt 2 



o \ 3 / 2 

2m \ 3/2 



(35) 



and in the high-field limit, Eqs. ^ and ( |l9| ) yield, 



,(3/2; 



(7) 



/2rr/ 



• rn 3,,2\ I/ 2 



(36) 



Both of these limiting cases are in full agreement with-jthc 
results expected for a 3D homogeneous electron gas.0 

The kinetic energy density is most easily obtained by 
using the second form of Eq. J22] ) with d = 3/2, viz., 



/ m \ 3/2 
~ \4n) 



Tj-l/2^5/2 [27^(3/2) _ Jp(5/2) 



(37) 



The functions T^/ 2 \^f) and T^^^ij) are readily ob- 
tained from Eqs. (|10|), and (pd|), and therefore 



4(2 7 ) 



5/2 



(27) 
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(38) 



The last four terms in Eq. ( p8[) are again quantum effects 
and are the origin of the de Haas- van Alphen effect. 

Using Eq. ( gjj ) in Eq. (|2^) yields the kinetic energy 
density in the zero-field limit 



£kir 



and in the high-field regime, Eq. ( 

Huj 27T 4 ?l 4 



£kin — — P ' 



(39) 

in Eq. ( Bgh gives 
(40) 



Equation (39) is the well known TF ijesult for a degen- 
erate noninteracting 3D electron gasja while Eq. (MG) 
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agrees with the result of Yonei et aZ.Jll and Banerjee et 
al.c5 Note that the functional form of the kinetic energy 
density is very different in the low and high-magnetic 
field limits. In particular, the p 3 term appearing in Eq. 
( [i"o| ) is a reflection of the one dimensional nature of the 
system along the field direction. 

In order to evaluate the exchange energy density, we 
require explicit expressions for the single-particle density 
matrix. In principle, Eq. ([32]) provides such an expres- 
sion. However, it is rather difficult to obtaip-the exchange 
energy for arbitrary magnetic field in 3D,li-3 so we focus 
here only on the the zero and high-field limits. 

In the zero-field limit, we set sinh(s) » s, and 



coth(s) 



Substituting these expressions into Eq. 



(|2|), and performing the inverse Laplace transform, we 
obtain 



P(r,r') 



1 



cosh(2y/aeJ) 
=sinh(2 y/aef) , 



(41) 



where a = {m/2irfi 2 ) | r — r' | 2 . This expression leads 
back to Eq. ( |35| ) upon taking the limit r — > r'. Equation 
(p0|) can now be used to yield the exchange energy func- 
tional. This is a straightforward, but tedious calculation, 
and since it only leads back to the well know expression 
for a. 3D free-electron gas, we simply quote the result 
here 



3 /3 



fox 



1/3 



e y /3 



(42) 



The high-field limit for the single-particle density ma- 
trix is obtained by replacing sinh(s) l/2exp(— s), and 
coth(s) w 1 in Eq. (0). After performing the inverse 
Laplace transform of the canonical density matrix, we 
obtain 



p(r,r') 



2hu 



m 



3/2 
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\/27r 2 \2wti A J I : - 

x sin(y / 2(e / - hu/2) \z-z'\) 

( .raw , , \ 
x exp^-i — (xy-y x)J 

x exp 



(-^x - x'f + {y - y') 2 )) ■ (43) 



(The high-field limits of the number density can be ob- 
tained from Eq. (43|) by taking the limit r — > r'). The 
exchange ener gy d ensity can now be calculated by sub- 
stituting Eq. (^) into Eq. ([so]) and transforming over 
to cylindrical coordinates: 



Ae 1 
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muj\ 2 



dzz 2 sin 2 



An^p 2 
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[°° drr(r 2 +z 2 )- 1 / 2 ex V (-^) 
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Fortunately, Danz and GlasserEl have already considered 
a similar calculation, and after several manipulations in 
accordance with their work, we obtain 



e /mu 

4^3 KIT J 

dx x\n(x 2 + l)exp(— pa; 2 ) 



p j dx a;ln(x )exp(— px ) 
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2p 



/>oo 
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Jo 



where 



8TT 4 h 



4*3 



-P 



(45) 



(46) 



The asymptotic behaviour of Eq. ( J45| ) in the low-density 
limit is given by 



(47) 



where C is Euler's constant. In the high-density regime 
(where the semiclassical TF theory is applicable), Eq. 
( B ) has the form 



_e*_\ /mw\ 2 i/2 

8^ 3 1 \ h J v 



(48) 



Similar to the kinetic energy density functional, the be- 
haviour of the exchange energy density is markedly dif- 
ferent in the low and high magnetic field limits. 



B. N=2 

Since d is integral, we can make use ofEq. (H) to ob- 
tain the electron number density in terms of elementary 
functions, viz., 



\2ttK) 



27 + E 



2(-l) fe . 



fc=i 



irk 



sin(27rfc7) 



(49) 



The low and high-field limits of the number density are 
given by 



m,Ej 
nh 2 



and 



P { '(l) = —r 
irh 



(50) 



(51) 



respectively. 

The kinetic energy density can also be expressed in 
terms of elementary functions, viz., 
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c ki] 



47T 



2 7 ^ (1) ( 7 )-^ (2) ( 7 ) 
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47T 

2(-l) fe 
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cos(27rfc7) 



2m 
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where we have used the fact that H^fa) = 1/2 V7. 
Therefore, in remarkable contrast to the 3D case, the 
the 2D kinetic energy density has exactly the same func- 
tional form regardless of the magnetic field strength; all 
of the de Haas-van Alphen type behaviour is buried in 
the 9^(7) dependence. This result is peculiar to 2D, 
since in all other dimensions, flw(7) is generally a func- 
tion of lj. For example, in four-dimensions, we have from 

Eq. m 



in the zero-field limit, and from Eq. (p7|) 



£kii 



ir 2 h 3 2 
m z uj 



-P , 



(53) 



(54) 



in the high- field limit. 

The e xchan ge energy is calculated along the same lines 
as Sec. [II A. The single-particle density matrix for ar- 
bitrary magnetic field is obtained by setting d = 1 in Eq. 
(p2n, and is given by 



p{r,v ) = — -exp 



2h 



-(x'y - y'x) ) x 



(- 

(i) frruj 



exp 
L 



raw 



r' I 2 l x 



(55) 



where we have summed only up to tha Fermi level 71/ , 
and used the finite summation relationtil 



e 2 /mui\ 3 / 2 
"n \~2h) 

T(n f + 



k=0 

2 + fe)r(2fe + l/2) 



T(n f + 1 - fe)r(2fe + 2)r(fc + 2)T(k + 1) 
F[-(n f - k), (2k + 1/2); (2k + 2); 2] , (58) 

where F(a, b; c; d) denotes the hypergeometric series.! 
The limiting forms of Eq. (pq ) in the zero and high-field 
limits can be made more transparent by first examining 
at the limiting forms of p(r, r') in Eq. (||). 
In the zero-field limit, we have 



P(r,r')= - 



1 



— kfj^kf I r-r' 



(59) 



where the Fermi wavevector kf = (2ms / /Ti 2 ) 1 / 2 . Once 
again, Eq. ( |50| ) is recovered by taking the limit r — > r' in 
Eq. (|59|). The exchan ge e nergy in this limit is given by 
an application of Eq. (pO|), viz., 



^2^/3 

3 r ' 



(60) 



This is the usual 2D free-electron gas result. In the high 
field limit, the single-particle density matrix takes the 
form 

muj f . mu> , \ 
9 (r, r = — ^exp -«r-r(a; y-yx)\ 
nn \ 2nn / 

(-£|((^--') 2 + (2/-2/') 2 )) , (61) 



x exp 1 



and upon using Eq. 



we find the simple result 

™2 

3/2 



2V2' 



(62) 



Equations (|60|) , and ([H]) should be contrasted with Eqs. 
([44), and (f45|) . Obviously, the key difference between 
the two sets of functionals is the universal form assumed 
by the exchange energy in 2D. Specifically, the exchange 
energy functional in 2D is explicitly independent of 7 in 
limit of strong magnetic fields. 



J2^ k Hx)=L^(x) 



n=0 



(56) 



The exchange energy in 2D can be evaluated from Eq. 
( ^p| ) , provided that we make use of the identity^ 



r(2 + n) ^ (2n-2fc)!(2fc)!4 2) (2a;) 

-\ t=^ u — Err - 



k=0 



T(2-k) (n-k)\ 



Upon substituting Eq. (55) into Eq. (p0[), we obtain 
after some tedious algebra 



IV. CURRENT-DENSITY FUNCTIONAL 
THEORY IN STRONG MAGNETIC FIELDS: 
THE TFD APPROXIMATION 

In this section, we take a pedestrian approach to the 
construction of a current density functional theory ap- 
plicable to inhomogeneous electronic systems in strong 
magnetic fields. The foundations of the. CDFT formal- 
ism have already been provided by VR£I, and we invite 
the interested reader to refer to this work for rigorous 
proofs and details. 
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Consider an inhomogeneous electronic system in a 
magnetic field. Within the TFD approximation, we as- 
sume that the non-uniform electronic system can be lo- 
cally approximated as a uniform electron gas character- 
ized by the density p{r), local Fermi wavevector fc/(r), 
and canonical (or "conjugate") current density j p (r) (see 
below). Thus, a local approximation to the energy func- 
tional is obtained by summing over all N-dimensional vol- 
ume elements, treating each N-dimensional volume ele- 
ment as independent; this is the local-density approxima- 
tion (LDA) in which the density is assumed to vary slowly 
on the scale of the magnetic length Ib — (frc/eB) 1 / 2 . 

The local approximation to the total ground state en- 
ergy functional for an N-dimensional many-electron sys- 
tem characterized by a scalar potential v(r) and vector 
potential A(r) is given by 



E T fu= / dr s[p(r), \ u(r) |] 



m 
~2 



dr 



Ml 

p{r) 



E n[p(r)} + dr p(r)v{r) 



de[p{r), | v 
dp((r) 



= M = Mki: 



(69) 



where the chemical potential \i now consists of kinetic 
plus exchange contributions to the free electron gas in 
the fictitious magnetic field. 

The minimization of Eq. (|6^ ) is done first with respect 
to j(r) (or equivalently j p , since p and A are held con- 
stant). This leads to an Euler-Lagrange equation for the 
equilibrium orbital current density, which may be cast 
in terms of the density p{r). A substitution of the re- 
sulting expression for the orbital current density into Eq. 
( |63| ) yields a functional, -Etfd, solely in terms of p[r). A 
minimization of this functional gives the Euler-Lagrange 
equation for the ground state density distribution of the 
many-electron system. 

Below, we treat separately the 3D and 2D systems 
within the CDFT, restricting our attention to the case 
of strong magnetic fields . Th e ener gy-density functionals 
obtained in Sees. |IIIA| and fllB now find relevance in 



(63) the construction of the total energy density e = ekir 



where En is the nonlocal Hartree contribution, and e — 
ckin + £cx is the total (gauge-invariant) energy density 
of the uniform electron gas in a fictitious magnetic field, 
defined by the vorticity v(r) 



i/(r) 



The physical orbital current density is given by 



where 



j(r)=j p (r) + — p(r)A(r) , 

mc 



j P ( r ) = — Im[Vp(r,r') | r=r /] 
m 



(64) 



(65) 



(66) 



The rationale behind considering the total energy e in 
a fictitious magnetic field is as follows: For a given lo- 
cal value of j p , we construct a fictitious vector potential 
A„(r) which produces j p in a uniform electron gas of den- 
sity p. The potential A„(r) couples only to the current 
and density. Since we roust have j = everywhere within 
a uniform electron gasjij we find 



A„(r) 



mc j p 
e p(r)' 



(67) 



The magnetic field associated with A„ (r) is nothing more 
than the vorticity defined in Eq. (64). The relationship 



between the fictitious magnetic field and the true mag- 
netic field is given by 



u(r) 



p{r)' 



(68) 



Therefore, the functional e is simply constructed by con- 
sidering the Euler equation for a homogeneous electron 
gas in a fictitious magnetic field v{r), viz., 



A. N=3 

The total energy density e is obtained from Eq. (|6 
with 



h I v 



4^4 



2 2 

-p -e 



8n 2 h 



(70) 



where the last term is the exchange contribution to the 
chemical potential. Using Eq. we immediately ob- 

tain 



e[p> I v I] = 2 I v I P 



27T 4 n 



4fc4 



/ m \ 
\8ir 2 h) 



3m 3 

1/2 , 



u | 1/2 p . 



(71) 



Notice that Eq. (^lHs formally equivalent to replacing 
w — *| v | in Eqs. (40() and ( p5|) (e.g., the total energy 
density in the fictitious magnetic field). As discussed 
above, a variation of the £tfd with respect to j gives 





de 


ivx 




m 


dv 


(-- 


A^h i p 2 


\2m 


TO 4 | V | 3 


e 2 / 


\ 1/2' 



2 V87T 2 ^ I V 



Vp x z 



In the limit of very strong fields, we obtain 



- — Vp x z 

2m 



(72) 



(73) 



which is exactly the result of Vignale and Skudlarski.Eil 
The vorticity can now be developed by using Eq. (J7S 



7 



in Eq. (p8|). Inserting the resulting expression for the 
vorticity into Eq. ( |7l|) gives all of the necessary com- 
ponents for the complete characterization of the ground 
state energy functional, Eq. (|63|). 

One can easily improve on the above TFD theory 
by including gradient— corrections (the so-called von 
Weizsacker correctionalll) in the kinetic energy density. 
The explicit inclusion of such gradient terms would take 
into account the magnetic-induced anisotropy of the 
density along, and perpendicular to the field direction, 
thereby giving a more realistic treatment of the system 
in a strong magnetic field. Such a calculation has al- 
ready been performed by Mazzolo et aL,c3 but without 
the exchange interaction. Their final result for the total 
energy functional is valid within the Thomas-Fermi-von 
Weizsacker approximation, and is given by 



E 



TFvW 



2 IVxH 2 



4*4 



8m 
8 m 



P 

4*4 

7r n 

m 4 uj 3 



2n 4 h 



37r 4tJ 
4 3 P 



V_Lp 



V±p 



v ll/° 



24m 



+ dr pv(r) 



(74) 



where Vj_ and V|| are the gradients perpendicular and 
parallel to the magnetic field direction respectively. Our 
formalism offers an improvement over Eq. ( |74] ) by in- 
cluding the effects of the exchange interaction. 



B. N=2 



e[p,\v\]=n { —)p~ 



2y/Y 



(77) 



Recall that in 3D, the magnetic-TFD theory is readily 
distinguished from the non-magnetic theory because the 
local 3D energy functionals enjoy very different behaviour 
in the low and high- field limits. Conversely, the local 2D 
energy-density functionals have no such signature. 

The current density is obtained with the aid of Eq. 
©, viz., 

1 ^ de 

J = - Vx 7T 

m av 








{ m 
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I 2nh 1 v 1 





Vp x z 



/_m_y 



Vp x z 



(78) 



In the limit of very strong fields we again recover Eq. 
([73|). MacDonald and Girvinllj have shown some time 
ago that Eq. ( |73|) is a special case of a powerful identity 
expressing the single-particle density matrix in terms of 
the density. Indeed, for a 2D electron gas in a uniform 
magnetic field, Eq. (|73|) is the exact orbital current den- 
sity distribution when only the lowest Landau level is 
occupied. 

The recent interest in the optical, and electronic prop- 
erties of the 2D electron-gas in a strong magnetic field 
(e.g., quantum-hall effected, 2D Wignec-crystallizationEil, 
edge reconstruction in quantum dotsE 2 ]) provides ade- 
quate motivation for the explicit calculation of Eq. ( [63] ) 
in two-dimensions. 

The vorticity is determined from Eq. ( f78| ) and Eq. 
(El, viz., 



The determination of the total energy density e in 2D 
is formally the same as in the 3D case. However, the 
simple substitution u> — >| v | in Eqs. ([32]) and ( |6l| ) for 
the present scenario is quite ambiguous. Obviously, this 
is because the energy-density functionals for N = 2 are 
explicitly independent of the magnetic field in the strong- 
field regime. So, unlike the 3D case, we have to be a 
little more careful when constructing the 2D total energy 
density e. 

First, we define the chemical potential for the homo- 
geneous 2D electron gas in the fictitious magnetic field: 



h I v 



7re / m 



1/2 



a/2 



2 2V2^2nh 
An application of Eq. (p3) yields 



U\v\ 7re 2 / m \ V2 , 



2y/2 \2irhJ 



(75) 



(76) 



If we now use Eq. ( |51| ) with uj — >| v |, Eq. (|76|) may be 

cast in terms of the density p, viz., 



u(r)\ 



i Tie f m \ 
2m A\f2m \2-kTiuj ) 
V^| 2 



1/2 



(79) 



where terms of the form V^p are ignored, as they eventu- 
ally integrate to zero. Using this expansion in Eq. fl76|), 
and with Eq. (|7Sl) for the current, we obtain 



m I j 



huj ire 2 ( m \ 2 , /9 



= - p 



H 2 e 2 



8m 

V^p| 2 



(—)' 

8 \2muj J 



7re 



P 



(80) 



Using Eq. (|80j) in Eq. (63), we finally obtain for the total 
energy functional, -Etfd 



3 



TFD 



dr 



(- 



2V2 



^l 2 p 



8 



h 2 e 2 / irh y ^ 
8m ~ ~8 \2muo ) + Mfiuo 

v(r) P (r) + - I dr'^^ 



(81) 



Notice that the first two terms have the same form as 
the usual zero-field TFD energy functional (the kinetic 
term is exactly the same as the zero- field TF term). Ex- 
tending Eq. d8l| ) to include gradient corrections in the 
kinetic energy is not possible because the, analogous 2D 
von Weizsacker term vanishes identically^ It is interest- 
ing to note, however, that the current density naturally 
provides a von Weizsacker-like (vW) gradient correction, 



viz., 



2m 



dr 



4p 



(82) 



with a coefficient that is l/|9-th the magnitude of the 3D, 
zero-field TFDvW theory.0 It is well known that the von 
Weizsacker term provides the smooth exponential fall-off 
for the density distribution as r — * oo. 

Equation ( pl| ) is applicable to any smooth external 
confining potential v(r). For example, choosing v(r) = 
l/2mwor 2 in Eq. ( pl| ) would provide a scheme for calcu- 
lating the ground state density distribution of a parabol- 
ically confined quantum dot in a strong magnetic field. 



V. CONCLUSIONS 

We have presented a general approach for the calcu- 
lation of the exact local-density functionals for an N- 
dimensional electronic system in a magnetic field. We 
have found that 2D systems enjoy a preferred status over 
other dimensions, in that the form of the 2D kinetic and 
exchange energy functionals are explicitly independent 
of the magnetic field. We have also provided an outline 
for the construction of a CDFT for 3D and 2D systems, 
which illustrates how the local-density functionals ob- 
tained in this work may be used to describe non-uniform 
electronic systems in a strong magnetic field. The results 
presented here may be of use in a wide variety of problems 
of current physical interest. In particular, our results for 
2D systems should be of relevance to researchers inter- 
ested in density-functional methods for low-dimensional 
electronic systems. 

Although the present work has focused on spin- 
averaged densities at zero temperature, an extension to 
include spin and finite temperatures is straightforward. 
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